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We investigate the effects of possible decay of neutrinos from a neutrino factory in a long 
baseline experiment. We consider the neutrinos from a factory at CERN and the detector 
to be the 50 kTon iron calorimeter (ICAL) detector proposed for India-based Neutrino 
Observatory (INO). We found considerable depletion of muon yield at INO for certain 
value of decay parameters. 



High energy neutrino beams are now playing a very crucial role to determine the 
yet unknown features of neutrino physics. Apart from determination of neutrino 
mass hierarchy and determination of 8 13 based on neutrino oscillation theory, another 
possibility could be to search neutrino decay over long baselines. The neutrino decay 
scenario has been discussed earlier in Ref. [H El El Hj . The neutrino decay scenario 
has also been discussed in connection to atmospheric and solar neutrinos in Refs. 
El [TJ, [HI El [TOl [HI [12]. Decay of Supernova neutrinos has also been adressed by 
several authors [131 HS1 ESI HZ]- The possible neutrino decay for cosmic ultra 
high energy neutrinos are also studied [18] and also neutrino decay is addressed in 
the context of recently proposed unparticle in [TQl [2U1 ETJ . Furthermore, MiniBooNE 
collaboration [22] searching for — > v e oscillation recently reported null result within 
the energy range 475 < E u < 1250 MeV and hence disfavours two neutrino oscillation 
theory, however, excess of v e events has been observed within the energy interval 300 < 
E v < 475 MeV. Attempts have been made to explain such result either through the 
inclusion of two sterile neutrinos [23] or through neutrino decaying into unparticles, 

eh. 



In the present work, we explore the possibility to observe neutrino decay for neu- 
trinos from neutrino factory through longbaseline neutrino experiment assuming a 
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beam from muon storage ring at CERN to India-based Neutrino observatory (INO) 
covering the distance of around 7152 Km. The detector at INO site is a magnetized 
Iron Calorimetre (ICAL) [21] which has an unique ability to distinguish charge of the 
particle passing through it. The mass of the proposed is around 50 kTon An exposure 
time of five years will be enough to collect large number of events. Assuming such 
mass and exposure time, in the present work we studied the 

Neutrino beams are generated from the decay of muons, /i ± — > e ± + v e (u e ) + u ^{v^) 
from the srtaight section of the muon storage ring [25l [26]. For unpolarized muon 
beam, the flux distributions of v e and are given by 




itL 2 El v » \ E u 
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El (E, - 2E Pe ) (2) 



where gi a b is the number of muons produced, J is the Jacobian factor arising due to 
transformation from rest frame to lab frame and is given by 

J = 7 (1 -0 cos 0) (3) 

E^ be the muon energy and E„ and Ep e are energies of produced corresponding 
neutrinos, 7 is the boost factor and j3 = p^/ E^ and 9 is the off axis angle which we set 
zero. For our analysis, we set all those parameters as gi a b = 0.35 x 10 20 considering 35% 
efficiency of the produced muon number with energy = 20 GeV. The parameter L 
is the length traversed by the neutrino from the source to the detector through earth. 
In the present work, we take L = 7152Km which is the distance between the source 
at CERN to the detector at INO site at PUSHEP (11°5'N, 76°6'E). The distribution 
of flux for both and v e are shown in Fig. 1 and Fig. 2 respectively. E v varies 
from 0-15 GeV E Pe varies from 0-10 GeV and and we observe for both the cases a 
large number of muons hitting the detector. 

Next, we consider the path of neutrino traversing through the earth from CERN 
to the proposed INO site at PUSHEP considering earth matter density profile and 
assuming three flavor neutrino oscillation without CP violating phase. The CERN- 
INO baseline length of 7152 Km is very close to the "magic baseline" that produces 
null CP violation effect. 
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We consider v e — > as well as v e — > mode including three flavor oscillation. The 
cross-sections due to interaction of v with the ICAL (Iron Calorimeter) detector pro- 
posed at INO site are mainly arisinig through quasi-ealstic (QE) and Deep inealstic 
scattering (DIS). At low energy we have considered QE mode (E v ~ lGeV), and 
for E v > lGeV dominant contribution will come from DIS. We have also included 
resolution function of the detector obtained from exact simulation of the ICAL de- 
tector which correlate the energy of the incident neutrino on the the detector to the 
produced muons inside the detector [21]. 

In the present decay scenario, we consider that the neutrino states ^2) and \us) 
are unstable and they decay into the stable state of \ The exponential decay term 
for the ithe neutrino state is given by exp(— AitL/X^) where L is the baseline length 
in kilometers and 

\ di = 2.5km-^-— (4) 
CeV aj 

where oti = wij/r, being the mass of the neutrino mass eigenstate and r being 
the decay lifetime. 

The purpose of this work is to investigate whether the effect of any possible decay 
of neutrinos can be detected by an iron calorimeter detector sue as INO. For this 
purpose we have chosen the long baseline neutrinos from a future neutrino factory at 
CERN. 

We have considered the oscillation of such neutrinos with decay. The neutrino 
flavour and mass eigenstates are conencted by the MNS mixing matrix as \u a ) = 
Uai\vi), where U a i are the 3x3 MNS mixing matrix of the flavour and mass eigenstates 
of neutrinos. Here, a = e, /i, r (the flavour indices) and % = 1, 2, 3, the mass indices. 
In order to incorporate the effect of any possible decay of the mass eigenstates \v<^ 
and 1 1/3), we fold the matrix elements U a i with exponential decay terms such that 
Uai -> U ai exp(-4:nL/\ d J for i = 2,3. 

The matter effect induced due to the passage of neutrinos through earth matter 
for traversing the CERN-INO baseline is also taken into account. We have taken an 
average matter density of 4.14 gm/cc for the purpose. This average density for the 
particular baseline considered is calculated using PREM [27] earth matter profile. 

We have made a three flavour calculation and the oscillation parameters used in 
the calculations are 

Am 2 32 = \m\-m\\ = 2.21 x l(T 3 eV 2 , Am^ = \m\ - m\\ = 8.1 x 10" 5 eV 2 
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for two mass square differences and 



0. 
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45°, 12 = 33.21°, 9 
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9°. 



for three mixing angles. 

We calculate the variation of total muon yield at INO with the decay constant a 
for such a decay scenario, in two cases; a) variation with a 2 for different fixed values 
of a?, and b) variation with a 3 values for different fixed values of a 2 . The results are 
shown in Figs 3 and 4 respectively. From Fig. 3, one sees that decay effect due to 
a 2 is most effective in the region 0.0001 ^ a 2 ^ 0.001. From Fig. 3, it is also seen 
that the total muon yield depletes to even less than half as a 3 varies from 10 -6 eV 2 
to 10~ 2 eV 2 . With 0:3 = m 3 /r eV 2 and m 2 ~ \J Am 2 2 ~ 0.05 eV, the lifetime r can 
be estimated as ~ 10~ 14 sec. for a 3 = 10 -2 eV 2 and ~ 10~ 10 sec. for 03 = 10~ 6 eV 2 . 
This is consistent with the limits given as in Ref. [6]. Fig. 4 shows similar plots for 
variation of total muon yield with a 3 for three different fixed values of a 2 . From Figs 
2 and 3 it is seen that the decay effect is more sensitive to a 3 than a 2 . Some of the 
muon yield values are tabulated below for reference. 



OL 2 


a 3 


Total Muon 


(eV 2 ) 


(eV 2 ) 


yield 


10~ 6 


10~ 6 


525326 


10~ 6 


10' 2 


227087 


10- 2 


10~ 6 


425453 



Table 1. Muon yield at INO from a neutrino factory beam at CERN for different values 
of neutrino decay parameters a 2 and 03. There can be more than 50% depletion in muon 
yield at INO due to neutrino decay. 

From the range of values of a 2 and a 3 that produces significant decay effects as 
demonstrated in this calculation, we make an attempt to verify whether the sum of 
the masses nil + m 2 + m 3 for three neutrinos is within the cosmological limit. In 
doing this we take the range of decay parameter a as 10~ 6 < a 2 ^ < 10~ 2 . We have 

m 2 /r _ m 2 
m 3 



Oil 



m 3 



m 2 



m 3 /r 



(5) 




Am| 2 



(6) 
(7) 
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mi = ^m\- ISm\ x (8) 
M = mi + m2 + rriz (9) 

Eqs. (5 - 9) is evaluated for normal hierarchy only. Using Eqs. (5 - 9) we compute 
M and found that within the ranges of a considered above M is always less than the 
cosmo logical bound of 0.7 eV. 

In summary, we have investigatd possible decay of light neutrinos in a long base- 
line experiment for a neutrino beam from CERN neutrino factory to ICAL detector 
proposed for INO covering a baseline length of 7152 Km. A unique facility of ICAL 
detector, the charge identification, has opened up the possibility of such study. A 
significant depletion of neutrino flux caused by the neutrino decay can be observed in 
the form of depleted muon signal at INO for a reasonable choice of decay parameters 
QJ2) «3 and other neutrino oscillation parameter inputs from atmospheric and solar 
neutrino experiments. 
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Figure Captions 



Fig 1. The flux from a neutrino factory vs different energies (E v ) for initial 
muon energy = 50 GeV. See text for details. 

Fig. 2 The z/ e flux from a neutrino factory vs different i7 e energies (E„ e ) for initial 
muon energy E^ = 50 GeV. See text for details. 

Fig. 3 The variation of total muon yield at INO with decay parameter a 2 for three 
different fixed values of decay parameter 0:3. 

Fig. 4 The variation of total muon yield at INO with decay parameter a 3 for three 
different fixed values of decay parameter a 2 . 
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